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1.0 INTRODUCTION MD SUMMARY 


Business and Technological Systems, Inc. has installed a program in the 
ATSOCC on-line computer which operates with attitude sensor data to produce a 
smoothed real-time orbit estimate. This estimate is obtained from a Kalman 
filter which enables the estimate to be maintained in the absence of T/M data. 

This report describes the results of analytical and numerical investigations 
into the sensitivity of Control Center output to the position errors resulting 
from the real-time estimation. 

The results of the numerical investigation, which used several segments 
of ATS-6 data gathered during the Sensor Data Acquisition run on August 19, 1974, 
show that the implemented system can achieve absolute position determination with 
an error of about 100 km. Implying pointing errors of less than 0.2° in latitude 
and longitude. This compares very favorably with ATS-6 specifications of approxi- 
mately 0.5° in latitude- longitude. 
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2.0 "SENSITIVITY ANALYSIS 


While S/C position is a very essential piece of information in computing 
S/C attitude or pointing vector, these output variables are relatively Insensitive 
to changes in S/C position. Primarily this occurs because of the 6.62 ratio be- 
tween the radii of the orbit and the Earth. In the current study, it was necessary 
to consider several distinct sensitivities: 

a) Sensitivity of subpoint to S/C position. 

b) Sensitivity of S/C attitude to S/C position. 

c) Sensitivity of pierce point to S/C position. 

d) Sensitivity of "sensed" S/C position to errors' in sensor data. 

e) Sensitivity of S/C pointing vector to S/C attitude. 

One of the purposes of the current task was to create a program for the POP— 11 
which would mechanize the triangulation for position-sensing and run a Kalman filter 
to smooth these "sensed" values. A second purpose- was- to perform a study evaluating 
the effect on S/C attitude and pointing vector of the position errors incurred by 
this real-time orbit determination. Clearly all of the sensitivities mentioned 
above are relevant to these purposes. 

For that reason it appears valuable to collect what analytic information is 
available on these sensitivities. 

For these calculations we take the radius of the Earth to be 

R'~ 6378 km 

and the nominal ATS-6 orbit radius to be 

r = 42164 km. 

2.1 .Sensitivity of 'Subpoint to S/C Position 

Except for Earth flattening effects,, the sensitivity- of subpoint to position 
is given by the ratio of radii. That is, sensitivity to changes along the radius 
vector are zero; sensitivity of subpoint position to changes normal to the radius 
vector is 1/6,62 = 0.151. 
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( 2 . 1 ) 


3x 

SC 

The relations between position and central angle are 

3x 

= 111 km/deg (2.2) 


dx 

= 736 Wdeg . (2.3) 

2.2 Sensitivity of S/C Attitude and Pierce Point to S/C Position 

The difficulty with this question is that it is not well posed. There are at 
least three modes in which the S/C position could be changed. The S/C attitude 
could remain constant in inertial space; it could remain constant in the (changed) 
ARS; or the sensor output could remain constant. All of these are analyzed here; 
however, it appears that, in general, constant inertial attitude is irrelevant; 
constant ARS attitude would he the applicable analysis if the computer (or S/C) 
were operating in an open- loop mode with no sensor data; constant sensor output is 
the correct model for normal operation. 


(1) Constant Inertial Attitude 

The S/C is at (x,0,0) and pointing at the Earth center. An increment in y 
is given but the S/C attitude is unchanged. The pointing vector now intersects the 
Earth at a new point displaced by the angle 0, where 


We thus have that 


d0_ ^ 1 

dy R cos 6 

Thus, the sensitivity of pierce point latitude or longitude at the sensitivity of 
S/C attitude -with respect to the ARS to cross track or downtrack position errors 


*ARS denotes the Attitude Reference System, which for ONATT is an orthogonal, 
S/C-centered coordinate system with z pointed to Earth center, y pointed south, 
and X east. 
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is' approximately 


0.16E-3 rad /km 


or 


0.009® (lat/long)/km . (2i4) 

(2) Constant Attitude t-JR! ARS 

The S/C is at (x,0,0) and pointing at the Earth center. An increment in y 
is given and the S/C rotates so as to remain pointing down the local vertical. 

The pointing vector now intersects the Earth at*a new point displaced by the angle 
0, where 


tan 0 = . 

X 


We thus have that 

2 

d cos 0 

dy X * 

Thus, the sensitivity of pierce point lat/long with respect to cross track or 
downtrack position errors is approximately 

4^ = 2.3 E-5 rad/km 
dy 

or (2.5) 


1.36 E-3" (lat/long) /km . 


(3) Sensor Output Constant 

When this mode is considered, some of the difficulties of the problem become 
very clear. Suppose that Polaris lay on the polar axis and the ESA/PSA combination 
were used for, attitude determination-. Then changes in S/C position, together with 
attitude changes keeping sensor output constant, would produce no change in attitude 
with respect to the ARS-. All changes in the pointing vector, then, would occur 
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Fig. 2 

Case (2) Geometry 
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Fig. 3 - 

Case (3) Geometry 
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from changes in the subpoint and (2) and (3) are identical for such an ESA/PSA. 

In actual practice, Polaris lies off the polar axis and therefore a change in 
S/C position can change the attitude in yaw. Through coupling, this can then 
change roll and pitch attitude. These are second order effects, however, and when 
ESA/PSA is used for attitude control, it seems proper to say that the sensitivity 
of attitude to position change is zero and that sensitivity of pierce point to 
position is as given in 2.1. 

Consider the S/C at (x,0,0) pointing at EC and an illuminator at longitude 
a. An increment in y is given and the S/C rotates so as to keep invariant the 
angle to the illuminator. The pointing vector now intersects the Earth at a new 
point displaced by the angle 0. 

The global relationship between y and 0' is exceedingly complicated and since 
we are only interested in local sensitivity, let us deal with incremental values. 
Using the law of sines and an intermediate angle 0, we find that 

sin 6 

Differentiating this and evaluating at y 



The differential change in central angle, Yj is of course. 


dy = — dy . 

X 


The angle, w, between the LOS to EC and to the illuminator satisfies 


$ + n = 90 - Y > 


and to keep the sensor output constant, we require that 


n + p = 4* , 
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where p is the pitch angle. 

1 

It follows that 


P + (J)-p = 90-Y 


or 


dp = dp + dy 


Using the formulae now available, 

dp = 

Taking nominal values (for Rosman) 



of 

•j) = 3.771° 
rQ = 36329 km 


we obtain 


= -3.75E-6 rad/km 
= -2..15E-4°/km . 


( 2 . 6 ) 


This gives the sensitivity of attitude to position when using an interferometer 
at Rosman to perform attitude control. 

Defining the new pierce point by 

(R cos 0 , R sin 6) , 

we find that 

R cos 0 = X - £ cos (90-B-(j>) 

R sin 0 = y - £ sin (90-B“4>) 
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where i is the unknown distance from S/C to pierce point. These can be solved to 
obtain 


tan (90“3-<J) 


y~R sin 9 
x-E. cos 0 


Differentiating this gives 

• -dg.- -dy-Rcos0d0 y-Rsin6 
S ^ o sxn 6d0 . 

cos (90-3-(j)) ^ (x-Rcos0)^ 

The previously derived expression for d8 can be used and the formula evaluated at 
y = 0» 6 = 0, 6 = 90 - <J), to obtain 


d9 _ 1- _ (x-R)cos(^ 
dy ~ R Rr^ 

Using the numerical vlues previously derived for Rosman, we have 

~ = 1.568E-4 - 1.541E-4 
dy 

= 2.7E-6 rad/km (2.7) 

= 1.54E-4®/km 

Thus we see that when the interferometer is used for attitude control, there 
will be a change in attitude caused by position, but the sensitivity of pierce 
point is very small since this almost compensates the change in subpoint. 

These examples serve to illustrate the problems involved in giving a concise 
answer to the question of pierce-point or attitude sensitivity to position. Because 
DOC usually controls with ESA/PSA, this combination was used to generate the numerical 
studies described in Section 3. 

2.3 Sensitivity of "Sensed" S/C Position to Sensor Errors 

In theory, the sensitivity of sensed S/C positions to errors in sensor output . 
can be as much as 

■— = "10000 km/deg (2.8) 
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The sensitivity given here varies of course with S/C attitude, which sensor is in 
error, etc. However, the number is approximately correct. Some analytic examples 
VTxll be given below. 

Since the quantization of interferometer readings is 0.0014^ and that of ESA 
about 0.005®, the errors on a quantization basis alone range as high as 50 km. 
There is a serious problem with word length in the OCPS calculation, however, and 
these accuracies are barely achieved (see Section 3) either because of large, 
variations in sensor readings, or because of truncation in the PDP-11. 

A 

Example 1; The following analysis gives a figure of 9600 km/deg. Looking at 

Figure 4, suppose that the ESA is giving correct readings but the IF is in error. 

In order to preserve the ESA data and correct the IF data, the S/C estimate must 

move oh the line pointing S/C and Earth center. To determine ,the sensitivity, we 
dr 

must compute From the law of sines 

R sin Y = r sin a 

and therefore 


R sin (a+0) = r sin a 


From this we have that 


dr rcosg-Rcos (a+g) 

da sina 

Letting 3 = 22® as before, we have that a = 3.771® and 

~ = -9640 km/deg . (2.9) 

da 

Example 2: Another verification of the sensitivity arises from the following 

special case analysis: 

The approximate coordinates of Rosman, Mojave and Earth center are 
2.12E6ft -1.70E7 1.20E7 lat 35®12' long 277®7' 

7.73E6ft -1.52E7 1.20E7 lat 35®19' long 243®7’ 

0 , 0 , 0 . 
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Let us take a bypothetical S/C in equatorial synchronous orbit over longitude 
260“7' (“99“53'). Then by a rotation of coordinates first about z, then about y, 
and rounding numbers slightly, this configuration can be represented as 


Rosman 

(0, 

.5E7, 

2S7) = (0,a,b) 

a^ + b^ = 4.25E14 

Mo j ave 

(0, 

-.5E7, 

2E7) = (0,-a,b) 

Earth Center 

(0, 

0, 

0) 

S/C 

(x. 

0, 

z) 


X = 1.13E8 

x^ + z^ H r^ = !• 901E16 

z = .79E8 

+ (z-b)^ = 1.625E16 
rj = 1.6275E16 

Now let us suppose that the LOS vectors to Rosman and Mojave are measured 
correctly so that cos a ^2 also correct i 

= 2 tan ^ ^ ) 

/x^+(z-b)^ 

(nominal value, ci ^2 ~ ^*^92°) 

and suppose that there is an error in the z-component of LOS to the Earth center, 

but that the y-component is correct (=0). Then the S/C will have to change its 

2 2 

X and z coordinates (subject to constant x + (z-b) ) so as to obtain the observed 
values 


U3 23’ 


(nominal value, 02 ^ - 7.720) while holding a ^2 constant. 


In this very special case we have 


2 2 

X + (z-b) constant 
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or 


xdx + (z-b)dz = 0 


( 2 . 10 ) 


Also 


2 

a 


+'b^ 


2rir3COsa23 




SO 


(ra -rj^cosa23)dr3 


= r^radcosa23 


Since 


2 2^2 
= X + z 


we have 


tadr^ = xdx + zdz . . 


Using (2.10) we obtain 


Tadra = bdz 


so that 


dz 

dcosa 


23 


^1^3 

b (Xa-r^cosct^a) 


Evaluating this at the unperturbed position, we have 


l.ObElOft. (2-. 11) 

dcosc23 

i 

We now need to determine the sensitivity of cosaaa to errors in the z component 
of ESA LOS. Up to this point, the analysis has been atti tude- independent . Now the 
attitude must enter. We know that 
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cosa^j = + *2"2 '•3“3 


, and 


m. 


- 7i ^ 


2 ”3 


Since approximately, = 1, we can say that 

3cos«23 _ 

3llij *3' " ”3 

This quantity is almost independent of attitude, and is approximately the sine of 
the difference in S/C angle between Mojave and Earth center, i.e. 


We conclude that 


Bcosa. 


23 


9m„ 


sxna 


23 


0.13 


0 £ £ 7500 km/deg . 

2.4' Sensitivity of S/C Pointing to Attitude 


(2.13) 


For a fixed subpoint, the pierce point changes with S/C attitude. At the 


subpoint the total angle change on the surface is about 
change, i.e. 


r-R 

R 


times the space angle 



(2.14) 


As the pointing vector moves to the edge of the disc, this sensitivity increases 
for two reasons. In longitude, the convergence of the meridians increases the 
sensitivity to pitch in the vicinity of the poles. For either latitude or longitude 
the sensitivity increases as the pointing vector approaches tangency to- the Earth. 
The- derivative of Earth central angle is given by 


d9 _ rcosa 
da Rcos(a+0) 


(2.15) 
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where a -is the S/C angle and 0 is the Earth central angle from the subpoint, 
a = 0 - 0 this gives 


da 


5.6 


and when the pierce point is at Rosman (a = 5.66°, 0 = 35^) then 


d0 ^ 
da 


^'3hen 


(2.16) 
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3.0 NUMERICAL RESULTS 


We have studied -carefully two data spans made during the ATS-6 sensor calibra- 
tion maneuvers on August 19, 1974. These covered the periods 

08 19 13:01 to 08 19 14:17 

during which time the pierce point went from S/C subpoint to Rosman to (0,-47°) 
with three interferometer calibrations and 

08 19 16:06 to 08 19 16:50 

during which time the pierce point went from (-15.6,-76.7) to subpoint to Mojave 
to (44,-138) to (-44,-138) with four calibrations. 

In Table 3 appear the pierce point coordinates at the relative times in the 
two intervals when interferometer calibrations were done. These calibrations are 
performed after the S/C has achieved equilibrium oscillation in low-jitter mode at 
the desired pierce point. They provide mark points allowing a rough calculation 
of attitude versus time. 

The coordinate system in which position is expressed is the GRS (Greenwich, 
Earth-fixed). Variations in ^^(z) are essentially cross-track, X 2 (y) radial, and 
x^(x) downrange. 
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Angle 

Error 

Time 

X 


<!> 

.078“ • 

14:00 

x^ = 460 

.480 at 13:50 

0 

.155“ 

13:28 

x^ = -846 

-846 at 13:28 


.018“ 

13:49 

"l = 

-846 at 13:28 

<i> 

.097“ 

16:06 

x^ = 560 

634 at 16:32 

e 

.166“ 

16:29 

x^ = -885 

-880 at 16:33 


.029“ 

16:30 

x^ = -885 

-880 at 16:33 


Table 1 

Maximum Euler Angle Errors 
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Angle 

Error 

Time 

Aa(x) 

Aa(6) 

Lat 

-1.17" 

1-4:01 

-0.62" 

-0.43 

Long 

-0.38° 

14:01 

-1.10" 

0.77 

Bear 

0.04" 

13:46 




Lat 

-1.06" 

16:25 

-0.77" 

-0.37 

Long 

0.46° 

16:39 

-1.16" 

0.90 

Bear 

0.06" 

16:40 




Table 2 


Total 

-1.05" 

-0.33" 

-1.14" 

0.26" 


Maximum Pierce Point Errors 



subpt cal at 750 sec*^ 

Rosman cal at 3000 sec. 
(0,-47) cal at 4000 sec. 
(-15.6,-76.7) cal at 200 sec. 
subpt cal at 850 sec. 

Mojave cal at 1600 sec. 
(44,-138 cal at 2100 sec. 


Table 3 


Pierce Points in Intervals 1 and 2 



3*1 Sensor Errors 


Certainly the first question which we wish to answer in this study is how 
good are the individual pseudo*-measurements of position. As we have seen in 
Section 2.3, the variations in sensed position will certainly be at least 50 km. 

On the other hand, any inaccuracies in sensor calibration, excessive noise in the 
sensor, or truncation in the PBP-11 will increase the variation. 

The plots of sensed positions versus disc ephemeris show several things: 

First there are the very large bias offsets. These biases are quite consistent 
in the two intervals, being about -750 km in and -450 km in The offset in 

^2 ranges between -40 and -100 km; however, there appears to be an attitude 
dependent component in all three components of about 100 km and in X 2 this is 
as large as the bias. * 

The second striking thing about the plots appears most readily in Figures 3 
and 6. This is the 200-300 km change in sensed values of x^ which take place when 
the interferometer is in the calibration mode. The calibration mode is a self- 
calibration capability in which the ends of the interferometer baselines are 
electronically switched. By comparing the unambiguous counts obtained in cal on 
with those obtained previously and subsequently in cal off, a partial calibration 
can 'be performed on-"line without interrupting data flow. If such compensation were 
performed, the bias in x^ could be reduced to about -250 km. A small effect of 
calibration can be seen in X 2 and this is also such as to reduce the bias offset. 
Unfortunately, no effect of calibration can be seen in the component with 
largest error. The effect of calibration compensation on x^ would be, in itself, 
sufficient to bring the latitude errors (Table 2) down to about 0.6^. 

In looking at the calibration bursts in Figures 3 and 6, two other phenomena 
are observed. There seems to be considerable attitude dependence appearing in both 
the offsets and the variance. In Figure 3 for instance, both the offset and the 
variance appear to be larger at Rosman than at subpoint. In Figure 6, the offset 
grows as the pierce point goes from the southern hemisphere through the subpoint to 
Mojave. In Figure 6, the scatter during the calibration at (44,-138) appears 
significantly larger than at the other points. 
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The questions which are raised here can be at, least partially answered by a 
thorough analysis of all the data taken during the sensor calibration run. Such 
an analysis could remove the obseirved biases and resolve whether the observed 
changes are attitude dependent or functions of time. 
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3.2 Filter Errors 


In Figures 7-18 are plotted the filter errors (x-x) and the filter residuals 

(x -x) for the two data intervals, 
s 

The primary purpose of examining the filter error is to evaluate its dynamic 
response and observe its variance about itself. Both features appear good in 
these runs. The time constant for the system is long (**12 min) but the response 
is quite smooth, with variations of less then 5 km. In an operational system 
there is no necessity for good step response so the constant gains used here 
(0.02, 0.0002) are probably close to the desired values. 

The component. Figures 7 and 10, shows classical underdamped step response. 
Notice that the data gaps disturb the response very little. 

The x^ response. Figures 9 and 12, is also quite satisfactory, showing good 
tracking of the time-varying values of sensor output. Notice the drastic effects 
of the calibration. 

The x^ component shows the attitude dependence of the readings very clearly. 
This component looks very irregular but this is caused by a combination of 
attitude- dependent effects, calibration errors, and the high resolution of the plot. 

Figures 13-18 are graphs of the filter residuals (x -±) which drive the filter. 

s 

The equilibrium value of this variable should be zero and its variance is a measure 
of the sensor noise. In all of these plots, it appears that the sensor noise has 
a standard deviation of about 50 km, with strong time or attitude-dependent 
fluctuations (Figures 14, 15, 17). 

The conclusion to be drawn from these graphs is that the ’^sensor” noise is 
approxin^tely as expected and appears to be reasonably well uncorrelated. There is 
also strong attitude dependency which careful analysis over the disc could probably 
remove. Neither of these are as serious as the large biases which can certainly be 
removed by calibration. 

Some statistical analysis of the residuals was done, but over the entire data 
interval which, because of the filter dynamic response gave misleading values of 
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the variance. When this analysis was performed- on the sensor errors, more realistic 
results x^ere obtained. These are presented in Table 4 and show that the statistics 
are quite consistent. The variance of is somewhat larger in interval 2 but this 
is partially because of the greater effect of calibration. If the calibration 
outliers had been removed then the standard deviation of 140 would have been 
reduced to about 110. 

.Table 4 also shows the time correlation of the signals. The paramter is 

the correlation coefficient of the signal with t and p 2 is the correlation of the 

2 ^ 

signal with the orthogonal component of t . Because of the long time interval 
these remove very little of the high frequency component. Instead they reflect 
time or attitude dependent trends in the signal. 
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Interval 

Variable 

Mean (km) 

a (km) 


• a (km) 

Pt2 

a^2 (km) 

1 


-749 

54 

.71 

38 

.16 

37 

1 


-50 

29 

.59 . 

23 

.16 

23 

1 

^s"^3 

-390 

82 

.30 

79 

.26 

76 

2 

Xs-Xi 

-723 

54 

.54 

46 

• .34 

42 

2 

^s“^2 

-50 

28 

.05 . 

28 

.32 

27 

2 

0) 

1 

-384 

140 

.36 

131 

.28 

125 


Table 4 


Sensor Noise Characteristics 
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3. 3 Euler Angle Errors 

Despite the extremely large values of filter errors, the Euler angle errors 
were usually quite small, as the analysis (Section 2.2(3)) predicted. The largest 
values observed in the two runs appear in Table 1, The times at which these values 
occur are very highly correlated with the times at which the maximum values of x 
occur. However, the Euler angle excursions in roll and pitch are only about one- 
seventh as large as the central angle errors implied by the position errors. 

Figures 19-24 show Euler angle errors (radians) for the two intervals. Notice the 
almost perfect correlation of 6^ with the filter error correlation 

of 66 with x^-x^. The yaw error, 6i|), seems to correlate very well with 
although the geometric reasons for this are not clear. 
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3 . 4 Pierce Point Errors 


In Table 2, a summary is presented of the largest pierce point errors appear-; 
ing in the two data intervals. Together with the lat-long errors have been listed 
the contributions to pierce point errors coming from subpoint shift and from Euler 
angle change. Note that in longitude, the errors caused by subpoint shift and by 
Euler angle error tend to cancel as occured in the example (Section 2.2(3)); while 
in latitude, the errors added. It is not clear whether this is caused by the S/C 
orientation, the station - S/C geometry, or the combination of position errors. 

This is however a rather important point in analyzing the expected values of pierce 
point errors over all possible sensor errors. For instance errors in 

had the same sign in both runs. Possibly the lat-long errors would change 
their behavior if these errors had opposite signs. 

These pierce point errors represent the final touchstone of on-line orbit 
determination accuracy. The mean value of errors of more than 500 km in estimation 
of and x^ are clearly too large; however, total absolute errors of approximately 
100 km (2 counts) appear to be attainable by proper sensor calibration, and these 
would bring' the pierce point errors down to less than 0.2“ which - corresponding 
to an 0.03“ S/C angle - is better than ATS-6 specifications. 

Future applications with more stringent requirements could be met also but 
might require more stability and better resolution in the attitude sensors. Certainly 
it appears that sensor noise is a simpler problem than proper calibration and 
assurance that calibration parameters are either constant or have well-defined 
dynamics . 

Figures 25-30 show pierce point errors (radians) for the two data intervals. 

There is an almost perfect correlation between the z component of filter error 

"S lat. 

There is also very good correlation between 6 long and 63 but neither of these 
show significant correlation with either the Euler angle or position errors. 
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4.0 CONCLUSIONS AND lUECOMMENDATIONS 


This study has demonstrated the feasibility of real-time orbit determina- 
tion from attitude sensors to an accuracy of about 100 fan (resulting in pierce 
point errors of less than 0.2°) for a S/C with three independent 2-axis sensors 
having error characteristics like the ESA and interferometer on ATS-6. 

Experience with ATS-6 real data indicates a requirement for sensor call- - 
bration in order to achieve these results. With current calibration constants 
there are extremely large position biases, on the order of 900 fan, leading to 
pierce point errors of about 1°. 

Sensor error (the difference between the position pseudo-measurement and 
GSFC definitive orbit) has three components: A high-frequency, uncorrelated 

noise'; a large bias, noted above; and a low frequency variation. These errors 
effect the estimation in distinct ways. The noise component ,. having a standard 
deviation in down-range and cross-range of about 60 fan (this is consistent with 
sensor variations of about 1 quantum) , was effectively removed by the sequential 
filter to a level of less than 10 km. The large biases can be easily removed 
by changing the ATTLMP calibration constants to a value determined by an overall, 
average, calibration. 

The low frequency variations, having an amplitude of about 100 km, are the 
errors which determine the absolute error figure given above. In the current 
system, which neither compensates for these variations nor solves for biases, 

•they are the factor limiting achievable accuracy. 

The most important change recommended is to the ATTLMP calibration constants, 
forcing them to give results more consistent with the definitive orbit. This 
would require no structural changes. 

The second change recommended is to modify the T/M reduction program so 
that the interferometer calibration mode is used to adjust the pertinent model 
parameters. This is a structural (programming) change. 

The evidence of this study shows that .these two changes will achieve system 
accuracy of about 100 km. To improve accuracy beyond this would require analysis 
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of data over the entire span obtained during Sensor Data Acquisition to evalu- 
ate attitude dependent effects and over an extended time period to evaluate 
temporal and temperature effects* 

To avoid these data reduction tasks, the feasibility of accurate real-time 
bias identification via an extension of the existing FILTER could be investi- 
gated. 

Recommended changes to FILTER in the present system are few. To guarantee 
all the accuracy possible in the triangulation, it would be desirable to operate 
in extended precision in the calculation of "sensed” position. More accurate 
orbit propagation models could be used, but seem to offer very little improve- 
ment in the context of 50 km errors. The system cannot be called operational 
until the initialization procedure is improved to remove initial transients. 

This problem is considered in detail in BTS-TR-75-18 describing the FILTER 
program. 
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